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Abstract. We present results of an extensive search for 
molecular gas in 14 high redshift (z > 2) radio galaxies. 
Radio galaxies with redshifts between 2.0 and 4.3 were ob- 
served during several sessions on two different single-dish 
telescopes (IRAM and JCMT) and two interferometric ar- 
rays (VLA and OVRO). No significant CO emission to a 
2a limit of a few times 10 10 Kkms -1 pc 2 was detected 
in any of these objects. The limits on the CO emission 
achieved indicate that, assuming a CO-H2 conversion fac- 
tor similar to the Galactic value, the mass of enriched 
molecular gas in these galaxies is less than 10 11 M Q . This 
suggests that distant radio galaxies may not be forming 
stars at the extremely high rates that have earlier been 
proposed, although they may still be forming stars at rates 
comparable to those of local starburst galaxies. 
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1. Introduction 

One of the key questions in the evolution of galaxies con- 
cerns the epoch of maximum star formation. Since in the 
Milky Way star formation takes place in molecular clouds, 
the presence of large amounts of molecular gas in distant 
galaxies would be an indication that they are undergoing 
prodigious star formation. For this reason it is important 
to find and study molecular gas in the most distant galax- 
ies. 

The first high redshift galaxy to be detected 
in molecular em ission lines was IRAS F10 214+4724 
at z = 2.286 ( Rowan- Robinson et al. 1991 ). Several 



CO lines were dete cted ( Brown fc Vandcn Bout 1991 
Solomon ct al. 1992a ) which would- indicate a molec 



ular gas mass of ~ ICFMqEI, dominating the 
dynamical mass. However, recently it was discov- 
ered that IRAS F10214+4724 is gravitationally lensed 
( Broadhurst fc Lehar 1995 , |Graham fe Liu 1995 ). Given 
the estimated lensing factor between 3 and 20, its true 
molecular gas mass may be only ~ 1O 1O M0. Since 
the discovery of CO emission from IRAS F10214+4724, 
confirmed molecular line emission has been detected 
in one other high redshift object, the also gravitation- 
ally lensed "Cloverleaf" quasar H 1413+117 at z = 2.5 



(Barvainis et al. 1994). In spite of the uncertain lensing 
factor (estimated ~ 7) this quasar must also contain a 
large amount of molecular gas. Noteworthy is that both 
the Cloverleaf quasar and IRAS F10214+4724 have an 
AGN-type spectrum (see Elston et al. 1994), i.e., a spec- 
trum with high ionization lines produced by a hard ultra- 
violet continuum. In spite of this ionizing continuum there 
is a large mass in molecular gas apparently shielded from 
the UV radiation so that the molecules can survive. 

In addition, there have been reports on possi- 
ble CO emission from damped Lya absorption sys- 
tems ( |Brown fc Vandcn Bout 1993| , [Frayer et al. 1994[ ), 
but more sensitive observerations could not confirm these 
(Braine et al. 1996, Wiklind and Combes 1994). 

The discovery of the strong CO emission from 
IRASF10214+4724 indicated that it might be possible to 
observe molecular gas out to the highest currently known 
redshifts and suggested that many high redshift objects 
might contain large amounts of molecular gas. High red- 
shift radio galaxies (HZRGs) usually have very bright Lya 
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1 Throughout this paper we assume a Hubble constant Ho = 
100 km s -1 Mpc -1 and a deceleration parameter qo = 0.5; val- 
ues quoted from other papers are, where necessary, tacitly con- 
verted to this cosmology. 
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emission, with integrated fluxes of ~ 10~ 15 erg s _1 cm~ 2 
(luminosities 10 ergs -1 ), and can extend over more than 
10" (~ 100 kpc), indicating the presence of a large reser- 
voir of gas in and around these galaxies. The importance 
of HZRGs for the formation and evolution of galaxies in 
general is shown by the fact that radio galaxies were al- 
most a thousand times more common at z ~ 2 than at 
the current epoch. The bright emission lines in the optical 
spectra of HZRGs are a factor of 5 — 10 narrower than 
those of quasars, facilitating an accurate determination of 
their rcdshifts, which is a prerequisite for deep mm- wave 
spectroscopy due to the limited spectral bandwidth of cur- 
rent instrumentation. 

It has been proposed that HZRGs are protogalaxies 
undergoing their first major burst of star formation (see 
Eales and Rawlings 1993, and references therein). At low 
redshifts, radio galaxies have molecular gas masses of the 
order of 1O 1O M0 ( Mazzarella et al. 1993). Such masses 
are comparable to those of known luminous infrared star- 
burst galaxies like Arp220 ( Solomon et al. 1990| ). High 
redshift radio galaxies have been suspected of having vig- 
orous starbursts associated with the propagation of the 



Table 1. Single-dish observing sessions 



radio jet (Rees 1989, Begclman & Cioffi 1989) as one of 



the explanations of the alignment of the optical, infrared 



and radio axes found in HZRGs (Chambers et al. 1987 



McCarthy et al. 1987), although alternatives have been 
proposed as well (cf., Sect. 4.2). The starburst model 
proposed by Chambers and Chariot (1990) requires the 
formation of 10 11 M Q of stars over 10 7 yr. Such star- 
bursts should be accompanied by the production of large 
amounts of heavy elements, giving rise to the formation of 
copious amounts of dust and molecules. Evidence for the 
presence of dust comes from polarization measurements of 
HZRGs (e.g., Tadhunter et al. 1992; Cimatti et al. 1993) 
and detections of continuum emission from dust at sub- 



millimetre wavelengths (Dunlop et al. 1994, [vison 1995) 



The detection of large amounts of molecular gas in such 
objects would be strong evidence in favour of the massive 
starburst hypotheses. 

In an extensive survey of ultrasteep spectrum radio 
sources we have found 29 radio galax- 



ics at z > 2 ( 


van Ojik et al.1994 




van Ojik et al. 1996a, 


Rottgering et al. 1995a, 


Rottgering et al. 1996), a signifi- 



cant fraction of the total number of radio galaxies known 
at such high redshifts. The large luminosities of the ob- 
served emission lines indicate that there is a power- 
ful source of ionizing radiation in such objects, albeit 
hidden from the observers. Any molecular gas would 
have to be shielded from this ionizing radiation (as in 
IRASF10214+4724 and the Cloverleaf quasar) to prevent 
the molecules from dissociating. The anisotropic ionizing 
radiation is most likely oriented in a cone along the ra- 



dio axis ( Antonucci 1993 ) . Thus most molecular gas might 
be situated outside this cone of ionizing radiation and is 
likely undisturbed by the radio jet, which may be respon- 
sible for the large velocity dispersions of the extended 



Session 


Date 


Telescope 


Time 
[hours] 


Comments 


1 


Nov 1993 


JCMT 


56 




2 


Nov 1993 


IRAM 


20 


variable conditions 


3 


Jan 1994 


IRAM 


28 




4 


Apr 1994 


JCMT 


24 




5 


Jun 1994 


JCMT 


40 


variable conditions 


6 


Aug 1994 


IRAM 


103 


variable conditions 


7 


Dec 1994 


JCMT 


88 




8 


Apr 1995 


IRAM 


55 


variable conditions 



optical emission line gas (Van Ojik et al. 1996b). Thus, 
we may expect any CO emission lines to be much nar- 
rower than the optical emission lines (which have typi- 
cally 1000 km s -1 FWHM) and similar in width to the CO 
line widths observed both in IRAS F10214+4724 and the 
Cloverleaf quasar (~ 300 km s -1 FWHM). We here report 
on our search for molecular gas in 14 HZRGs. 

2. Observations and reduction 

2.1. Transitions surveyed 

Our targets have redshifts ranging from 2.0 to 4.3. Taking 
into account the instrumentation available, this allowed 
us to survey CO transitions from J = 3^2 (345.8 GHz 
rest frequency) for the objects with z ~ 2 observed at 
IRAM or the Owens Valley Radio Observatory millimeter 
interferometer (OVRO) in the 3 mm window, to J = 9— >8 
(1036.9 GHz rest frequency) for the objects with the high- 
est redshifts observing at the JCMT in the 1.3 mm win- 
dow. With the VLA the CO J = 1— »0 transition can be 
observed at 25 GHz for the sources with the highest red- 
shifts. 

Strategies for detecting CO at high z have been 
discussed by Van der Werf and Israel (1996), who 
show that if conditions in these objects are similar 
to those in local starburst galaxies, the J=3^2 to 
J=6^5 lines are the transitions of choice in detec- 
tion experiments. Significantly, the line ratios found 
in IRASF10214+4724 and the Cloverleaf quasar are 
indeed similar to those in local starburst galaxies 
( |5olomon et al. 19924 Parvainis et al. 19941). 

With the exception of B2 0902+34 (Lilly 1988) 
and 8C 1435+635 (Lacy et al. 1994), which were 
i ncluded because o f their the rmal dust emission 
( |Chini fc Kriigel 1994j |lvison 1995|) , the sources observed 
were selected from the Leiden sample of HZRGs. One tar- 
get, 0211 — 122, has an optical spectrum similar to that of 
IRASF10214+4724, with strongly absorbed Ly a indicat- 
ing the presence of dust (van Ojik et al.1994). Four other 
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Table 2. Log of the single-dish CO observations, with the 
r.m.s. noise a in the spectrum given at a velocity resolution of 
100 kms" 1 . 
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targets were selected because they showed evidence for as- 
sociated Hi absorption (Van Ojik et al. 1996b; Rottgcring 
et al. 1995a) and the remaining targets were selected be- 
cause of their large Lya halos and high redshifts. Thus, the 
14 objects we observed are representative of our general 
sample of HZRGs. 



2.2. Single- dish observations 

We observed the sources in 8 different observing sessions, 
summarized in Tables 1 and 2. 



Four sessions were carried out with the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kca, Hawaii. We 
used the SIS receiver A2 tuned so that the redshifted fre- 
quency of the most suitable CO transition for each object, 
at 220 to 250 GHz, was in the upper sideband. The tele- 
-scope half-power beam width (HPBW) at these frequen- 
cies is about 20". As backend we used the 2048 channel 
.Dutch Autocorrelation Spectrometer with a bandwidth of 
750 MHz. The single sideband (SSB) system temperatures 
in good atmospheric conditions were 300 to 450 K. We 
chopped the subreflector at 2 Hz with a throw of ~ 60" 
and applied beam-switching. 

Four sessions were carried out with the IRAM 30 m 
telescope at Pico Veleta, Spain. Usually, two or three re- 
ceivers were used simultaneously, so that up to three CO 
transitions could be observed with high net observing 
efficiency. Two 512 channel, 512 MHz wide filter banks 
were used as backend for the 3 mm and 2 mm receivers 
and an autocorrelator backend for the 1.3 mm receiver. 
The HPBW was 25" at 3 mm, 17" at 2 mm and 11" at 
1.3 mm. The subreflector was chopped at 0.5 Hz with a 
throw of 120". System temperatures at IRAM under rea- 
sonable weather conditions were 200 — 400 K at 100 GHz, 
300 - 500 K at 150 GHz and - 1000 K at 230 GHz (SSB). 

The instantaneous velocity coverage of a spec- 
trum ranged from 500 km s" 1 (JCMT at 230 GHz) to 
1200 kms" 1 (IRAM 30m at 100 GHz). The CO lines from 
our HZRGs are expected to be faint and broad, probably 
of the order of 300 kms -1 , covering a large portion of an 
instantaneous spectrum. Thus, broad low level ripples in 
the baseline of the spectra are a serious limitation in de- 
tecting these lines. In order to minimize baseline problems, 
integrations were carried out at several different local os- 
cillator (LO) settings. By extending the covered baseline, 
residual baseline undulations should be better distinguish- 
able and not be mistaken for true CO emission. Usually 
the LO settings were at 0, +150 and — 150 kms -1 with re- 
spect to the calculated redshifted CO sky frequency, but 
occasionally we extended the baseline to 600 kms -1 from 
the nominal sky frequency. In this way we achieved a total 
velocity range of - 1200 - 2000 kms" 1 . 

Reduction of the spectra from IRAM was done with 
the CLASS (Continuum and Line Analysis Single-dish 
Software) reduction package, while the JCMT spectra 
were processed with the SPECX spectral line reduction 
package. The observations of each observing shift and each 
velocity setting were examined separately before summing 
them into a final spectrum. Scans with instrumental prob- 
lems (e.g., curved baselines) or taken in poor conditions 
(atmospheric opacity r > 0.2) were discarded. The fi- 
nal spectra were smoothed with a 100 km s -1 FWHM 
Gaussian (i.e., a factor of 2 — 3 smaller than the ex- 
pected width of any CO emission line), and resamplcd. 
The antenna temperatures were converted to Raylcigh- 
Jeans main beam brightness temperature outside the at- 
mosphere (T m b [K]). The corresponding flux density scale 
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Table 3. Log of interferometric CO observations, with the 
r.m.s. noise a in the spectrum given at a velocity resolution of 
100 kms^ 1 . 



Source 


z 


Telescope 


Line 


tint 

[hrs] 


a 

[mjy beam -1 ] 


1243+036 


3.5856 


OVRO 


4-3 


32 


1.6 


1243+036 


3.5856 


VLA 


1-0 


18 


0.6 


1435+635 


4.2576 


VLA 


1-0 


10 


0.3 



for a point source with the JCMT is 18.9 JyK 1 and 
with the IRAM telescope 4.7 JyK -1 at 3 and 2mm and 
4.8JyK+ 1 at 1.3 mm. 

2.3. Interferometric observations 

Several systematic effects which limit deep spectroscopy 
using single-dish radio telescopes are reduced or absent for 
interferometric measurements, since interferometers only 
record correlated signals and are thus not sensitive to arte- 
facts arising in separate elements. Two targets were there- 
fore observed with interferometric arrays to verify the real- 
ity of apparent emission features found during single-dish 
measurements. 

The VLA was employed to search for CO emission from 
two of the highest redshift objects, 1243+036 (z — 3.6) 
and 1435+635 (z = 4.3). The VLA has the advantage of 
a large collecting area and sensitive receivers, but the dis- 
advantage of a limited bandwidth (50 MHz maximum, or 
650 kms -1 at z = 4). Moreover, limitations in the IF sys- 
tem at the VLA constrain observations using a 50 MHz 
bandwidth to a discrete set of central frequencies, in se- 
quential steps of 20 MHz and 30 MHz. For our observa- 
tions we used a series of overlapping 50 MHz bands cen- 
tered around the nominal systemic redshift, with 8 spec- 
tral channels per band and two circular polarizations, syn- 
thesizing a total bandwidth of 130 MHz for 1435+635 and 
of 80 MHz for 1243+036. The observational cycle included 
a sequence of short (5 minute) scans at each frequency set- 
ting, with scans on a nearby phase calibrator before and 
after each source scan. 

The source 1435+635 was observed in the D-array 
(1 km maximum baseline) in May 1995. The phase cali- 
brator was 1435+638, and 3C 286 was used for gain and 
bandpass calibration. Four frequency settings were used 
implying a total redshift range searched for emission of 
4.2420 to 4.2732 (observing frequencies around 21.9 GHz). 
The source 1243+036 was observed during reconfiguration 
from C array (3 km maximum baseline) to D array in Oc- 
tober 1995. Only the D-array spacings were employed in 
the final analysis in order to limit problems with phase co- 
herence. The source 3C 273, which is located conveniently 
close to 1243+036 at an angular distance of only 5°, was 



used for phase and bandpass calibration, and 3C 286 was 
used for gain calibration. Three frequency settings were 
employed, implying a redshift range searched for emis- 
sion of 3.5784 to 3.5932. The observing frequencies were 
around 25.1 GHz for this source. This frequency is at the 
edge of the VLA K-band, leading to LO phase lock prob- 
lems on some antennas, and to a generally higher noise 
level (by a factor 2.5 or so relative to the K-band center). 
In the end, only 17 antennas were used in the analysis. The 
beam (FWHM) of the VLA (in hybrid configuration BnC 
for 1243+036 and DnC for 1435+635) was about 2.5" at 
22 GHz. This beam size is similar to the extent of optical 
continuum emission. 

Standard phase, amplitude, and bandpass calibration 
was employed using the Astronomical Image Processing 
System (AIPS) at the U.S. National Radio Astronomy 
Observatory, Socorro. Spectral channel images were con- 
structed for each source at each frequency setting, and 
visually inspected for emission. Also, spectra were ex- 
tracted at the position of the radio nucleus, and inspected 
for emission. In all spectra the theoretical noise level was 
achieved. One potential problem with high frequency ob- 
servations is phase coherence. To check the phase coher- 
ence, a continuum image of each source was synthesized. 
In both cases the continuum flux density observed was 
close to the expected value (determined by extrapolating 
a power law spectrum from lower frequencies to 22 GHz), 
and the spatial distribution agreed with that expected for 
each source. 

The radio galaxy 1243+036 was also observed at 
OVRO during December 1994 and January 1995 at 
100.5 GHz in rcdshifted CO(4-3). Four tracks were ob- 
tained in various configurations of the six telescopes, with 
interferometer baselines from 35 to 242 metres. During ev- 
ery track 1243+036 was observed for a total time of about 
8 hours. Phases were calibrated using 3C273 which was 
observed every 30 minutes. Once every hour, measure- 
ments of 3C273 were also done to correct the pointing. 
3C 273 was also used as a passband calibrator. The flux 
density scale was calibrated with measurements of 3C 84, 
3C 454.3, 0528+134, Neptune and Uranus. The backend 
was configured with four adjacent and slightly overlap- 
ping bands, each covering 128 MHz in 32 channels. The 
CO(4-3) line from 1243+036 was searched for in lower 
side band (LSB) at 100.5 GHz. System temperatures were 
typically 200 K (SSB). After standard editing and calibra- 
tion, AIPS was used to produce the final data cube at 
a spectral resolution of 32 MHz ^kms -1 ). The OVRO 
beam size at 100 GHz was about 2" 

The interferometric observations with the achieved 
noise levels are summarized in Table 3. 
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Fig. 1. CO J = 4-^3 (a) and J = 6-^5 (b) spectra of the radio galaxy 1243+036 observed with the IRAM 30 m telescope in 
January 1994. The spectra are centred at the redshift of Lye* from a low resolution optical spectrum, were smoothed with a 
100 km s -1 FWHM Gaussian and block-averaged into bins of ~ 85kms _1 . Both in the 3 mm (a) and the 2 mm (b) spectrum a 
feature was observed at an offset velocity of ~ 320kms _1 from the optical redshift. In spectra with different LO settings these 
features were present at a constant sky frequency as expected for an astronomical emission line. The same object was observed 
in April 1995 with the IRAM 30m telescope in CO J = 4^3 ((c) and (d)). The centre is now at a corrected redshift from a 
high resolution optical spectrum. We have extended the baseline by adding many overlapping spectra with LO settings up to 
600 km s -1 from the redshifted sky frequency. In spectrum (c) the data was block-averaged into bins of 65 km -1 . In spectrum 
(d), the data were first smoothed with a 100 kms -1 FWHM Gaussian and then rebinned. The feature seen in January 1994 
should have appeared at + 1050 km s -1 . No evidence for any emission feature is seen in (c) and (d), leading to the conclusion 
that the tentative result from the earlier observations (a) and (b) was spurious. 



3.1. Limits on CO emission from single dish observations 

A number of features were observed that significantly ex- 
ceeded the noise level in the rest of the spectrum and 
remained when the receivers were tuned to different off- 
set velocities. We first concluded that such features, at a 
level of a few mK, were real CO emission features. How- 
ever, none of these features could be reproduced in dif- 
ferent observing sessions. When re-observed, there was no 
evidence for any emission down to lower noise levels or 
new features appeared at a different frequency Examples 
of this situation are shown in Figs. 1 and 2. Thus the ob- 
served features in the spectra were not due to line emission 
but were instrumental artefacts. The ability to detect faint 
broad emission lines depends strongly on these instrumen- 



tal features in a spectrum and not only on the formal noise 
level. A discussion of these effects in the JCMT spectra 
has been given by Rottgering et al. (1995b). 

We determined the final r.m.s. noise level of the spec- 
tra at a velocity resolution of 100 km s -1 over the entire 
observed velocity range (thus also including any possible 
spurious features) and these noise levels are listed in Ta- 
ble 2. Because of the occurrence of spurious features, some 
spectra with relatively short integration times but no spu- 
rious features have lower effective noise levels than other 
spectra with longer total integration times but with spu- 
rious features. 

From the noise levels we estimate upper limits for 
the line flux from the objects. If we assume that a 
CO emission line has a Gaussian shape with a specific 
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Fig. 2. (a) Spectrum of 00(7-6) of the radio galaxy 0211-122 observed with the JCMT in November 1993. This spectrum 
includes only the data with atmospheric transparency r < 0.1, i.e. from 4 of the 6 shifts (where one shift is 8 hours). There 
appears to be a strong emission feature centred at — 150 kins -1 from the optical redshift. (b) Spectrum of 00(7-6) of 0211 — 122 
taken with the JCMT in December 1994, centred at the tentative emission feature from spectrum (a). We have extended the 
baseline by combining spectra with 5 different LO settings to better distinguish any emission features. The tentative emission 
feature is of (a) absent in (b). Both spectra have been smoothed with a 100 km s -1 FWHM Gaussian and block-averaged into 
bins of ~ 60kms -1 . 



FWHM, we estimate a 2a upper limit to the CO flux 
as 1.06 x FWHM x 2 x aioo where cri o is the r.m.s 
noise in the coadded spectrum at a velocity resolution 
of 100 kins' 1 . We assume a FWHM of 300 kms -1 as for 
the observed lines in IRAS F10214+4724 and the Clover- 



leaf quasar ( Solomon ct al. 1992a , Barvainis et al. 1994 ). 
From the line flux limits we derive a 2a upper limit on 
the CO luminosity based on the lowest CO transition ob- 
served under good atmospheric conditions, using the fol- 
lowing formula, which is independent of beam size (see 
Solomon et al. 1992a): 



L' co = 3.25 x 10 7 Sco&v ^ - b 2 s £>£ (1 + z) 



-3 



(1) 



where Sqo^ v is the integrated line flux limit in Jy kms -1 , 
^obs is the observing frequency in GHz, Z?l the luminosity 
distance to the object in Mpc and the line luminosity L' is 
expressed in Kkms -1 pc 2 . The limits derived for the CO 
transitions used are listed in Table 4. 

3.2. Limits on CO and dust emission from inter ferometric 
observations 

No emission was seen in any spectral channel for either 
of the two objects observed with the VLA and OVRO. 
We derived 2a upper limits for the CO luminosity for the 
interferometric observations in the same way as for the 
single dish observations. The limits are listed in Table 4. 

The VLA continuum maps derived from the observa- 
tions of 1243+036 (25 GHz) and 1435+635 (22 GHz) do 
show continuum emission. In 1243+036 there is a con- 
tinuum peak of 3.5 mjy beam -1 at the location of the 



brightest hot spot seen in 8 GHz radio imaging obser- 
vations ( van Ojik et al. 1996a ). The continuum emission 
of 1435+635 has a peak of 2.5 mjy beam -1 and shows 
structure similar to the 8 GHz radio image. The contin- 
uum flux densities are consistent with non-thermal emis- 
sion from the radio sources extrapolated using the spec- 
tral index determined between 8 GHz to 22 GHz. Thermal 
dust emission, as was recently detected from 1435+635 
( [vison 1995 ), would at a rest frequency of 115 GHz only 
be expected at the (iJy level (estimated from the spec- 
trum of Ivison 1995). The radio continuum is at the level 
of ~ 3 mjy at this frequency. Thus any emission from dust 
is negligible compared to the synchrotron radio emission. 

No continuum emission was observed in the 100 GHz 
continuum map derived from the OVRO observations of 
1243+036 to a 2a upper limit of 0.4mJy beam -1 . 



4. Discussion 

4-1. Limits on the molecular gas masses 

For Galactic giant molecular clouds the relation between 
CO luminosity and H2 mass has been determined as 



a = M(H 2 )/L' co = 2.2n^/T b = 4.78 



(2) 



( |5olomon et al. 1992a| ), where Tb is the intrinsic line 
brightness temperature (which is constant for thermal- 
ized, opaque lines from different rotational levels), njj 2 
is the H2 density in cm -3 , Af (H 2 ) is the H2 mass in M Q 
and L' co has been given in Eq. (1). The conversion factor 
a depends on the excitation conditions of the molecular 
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Table 4. 2a upper limits on CO luminosity and molecular gas 

mass 



Source 



Line 



line luminosity 
[Kkmt 



pc 



M(H 2 ) 
[10 11 M 



Limits from the single dish observations 



0211-122 


CO(3- 


-2) 


2.336 


< 2.4 x 10 10 


< 1.0 


0214+183 


CO(6- 


5) 


2.131 


< 4.1 x 10 10 


< 1.6 


0355-037 


CO(6- 


-5) 


2.153 


< 4.1 x 10 10 


< 1.7 


0448+091 


CO(6- 


5) 


2.040 


< 2.1 x 10 10 


< 0.8 


4C 41.17 


CO(4- 


3) 


3.800 


< 1.3 x 10 10 


< 0.5 


0748+134 


CO(3- 


2) 


2.410 


< 1.0 x 10 10 


< 0.4 


0828+193 


CO(3- 


2) 


2.572 


< 3.2 x 10 10 


< 1.3 


B2 0902+34 


CO(4- 


3) 


3.395 


< 1.3 x 10 10 


< 0.5 


0943-242 


CO(8- 


7) 


2.918 


< 2.6 x 10 10 


< 1.0 


1243+036 


CO(4- 


3) 


3.570 


< 1.3 x 10 10 


< 0.5 


1435+635 


CO(4- 


3) 


4.257 


< 1.1 x 10 10 


< 0.5 


1707+105 


CO(3- 


2) 


2.345 


< 1.8 x 10 10 


< 0.7 


2202+128 


CO(8- 


7) 


2.704 


< 2.3 x 10 10 


< 0.9 


2251-089 


CO(6- 


-5) 


2.000 


< 3.5 x 10 10 


< 1.4 



Limits from the interferometric observations 



1243+036 
1243+036 
1435+635 



CO (4-3) 
CO(1-0) 
CO(1-0) 



3.570 
3.570 
4.257 



< 1.0 x 10 



10 

< 6.5 x 10 10 

< 2.8 x 10 10 



< 0.4 

< 3.3 

< 1.3 



gas and may be either higher or lower in extremely lumi- 
nous and active galaxies such as IRAS F10214+4724 and 
HZRGs than in the Milky Way. More fundamentally, in 
order to be able to trace H2 with CO lines, the molecular 
gas must be enriched, and, while subsolar CO abundances 
are allowed, the CO column densities must be sufficient to 
produce optically thick 12 CO lines. Solomon et al. (1992a) 
show that for IRAS F10214+4724 the conversion factor a 
for the higher CO transitions is similar to the Galactic 
value of 4.78 for CO J = 1— »0. Using their value a w 4 for 
the lowest CO transitions observed under good conditions, 
we derive upper limits for the (enriched) H2 masses for our 
14 HZRGs of slightly less than 10 11 M . For consistency, 
we have used the same value for the CO-H2 conversion for 
the limits from the interferometric observations of CO(l- 
0), because it differs only little from the Galactic value 
of 4.78. The upper limits for the CO luminosity and the 
corresponding H 2 mass are listed in Table 4. 

The limits on the H2 masses that we find are 
marginally smaller than the molecular gas mass that was 
initially estimated for IRAS F10214+4724 (10 11 M Q be- 
fore it was realized that it is lenscd. In low rcdshift ra- 
dio galaxies, Mazzarella et al. (1993) found molecular gas 
masses of order 10 10 M Q . Thus, we find that radio galaxies 
at high rcdshifts are not much more gas-rich than those at 
low redshifts. As a prototypical gas-rich starburst galaxy 
and ultra-luminous IRAS galaxy (ULIRG) at low redshift, 



Arp 220 has a large molecular mass, M (H 2 ) ~2x 10 10 M 
(Solomon et al. 1990 and references therein). Although it 
is possible that our HZRGs contain amounts of molecu- 
lar gas similar to the ULIRGs, they cannot contain much 
more. 

4-2. Implications of these and other millimetre observa- 
tions at high redshift 

The two high z objects with confirmed detections of 
emission from molecular gas, IRAS F10214+4724 and the 
cloverleaf quasar, are both gravitationally lensed. This 
could explain why the present and other searches for 
molecular gas in high redshift galaxies and quasars have 
been unsuccessful so far and might be taken as a discour- 
agement for further searches. However, for the two ob- 
jects detected the lensing amplification is unlikely to be 
much more than a factor of 10-20 ( |Barvainis et al. 1994 ), 
which, although rendering the true CO flux density too 
low for current instrumentation, may be observable in 
the future. The amplification of the flux densities of 
IRASF10214+4724 and the Cloverleaf quasar imply that 
the underlying systems are not necessarily of extraordi- 
nary properties. Although still quite luminous and with 
~ 10 10 M of molecular gas quite gas-rich galaxies, they 
may well be "normal" active galaxies in the early universe, 
so that that many active galaxies, quasars and radio galax- 
ies at high z may contain similar amounts of molecular 
gas. 

The detection of submillimetre continuum emis- 
from a number of high redshift radio galaxies 



sion 



and quasars in several searches (McMahon et al. 1994 



Dunlop et al. 1994, Isaak et al. 1994, Ivison 1995) has 



been used to estimate gas masses using simple assump- 
tions on the dust temperature and gas/dust ratio. The 
estimates of molecular gas masses of 0.5 — 1 x 10 11 Mq 
from these continuum observations are not in contradic- 
tion with the non-detections of molecular gas in these 
objects and our sample of HZRGs. Far-infrared lumi- 
nosities of 10 13 — 10 14 Lq, as implied by the submil- 
limetre dust detections, imply star formation rates of 
10 3 — lO 4 M yr~ 1 if the FIR luminosity is powered by 
star formation. By extrapolating from the somewhat lower 
lumi nosity (Xfir > 1O 12 L ) ultra- luminous IRAS galax- 



ies ( Sanders et al. 1986 ), a value Lfir/A/h 2 of 100 — 
1OOOL /Mq is expected at these luminosities. The im- 
plied molecular masses are of the order of lO n M , only 
slightly more than the present upper limits. 

One of the explanations proposed for the strong align- 
mcnt of the optical, near-infrared and radio ax es of HZRGs 
( jChambers et al. 1987| , ^McCarthy et al. 1987| ) is that vig- 
orous star bursts are associated with the propag ation of the 
radio jet ( Rees 1989 , Begelman fc Cioffi 1989 ). In some of 
these radio jet-induced star formation scenarios, star for- 
mation rates of the order of lO 4 M yr _1 over a period 
of ~ 10 s yrs are invoked to explain the alignment effect 
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flChambers fc Chariot 1990| , phambcrs ct al. 1990| ). Such 
star formation rates would imply that a mass of order 
10 12 M© in (molecular) gas is converted into stars. The 
fact that we do not observe such quantities of (enriched) 
molecular gas in any high redshift radio galaxy is relevant 
to this scenario and can have two explanations: 



The HZRGs that we observe have already consumed 
most of their molecular gas and formed the bulk of 
their stellar population. Jet-induced star formation 
may have played a role, but we are observing the object 
in a later phase where the radio source has already ex- 
panded into a powerful Faranoff-Riley II (FR-II) type 
source with an aligned stellar component in the galaxy 
hosting the radio source. Mazzarella et al. (1993) found 
molecular gas in several low redshift radio galaxies de- 
tected with IRAS, but not in the most powerful FR-II 
radio galaxies. They suggested that the powerful FR-II 
galaxies might be in a later phase where most of the 
molecular gas has been converted into stars. 
The alignment effect may not be due to jet-induced 
star formation, but to another mechanism (cf., Mc- 
Carthy 1993 for a review) such as scattering of contin- 
uum radiat ion emitted along the radio axis from a hid- 
den AGN (|Tadhuntcr ct al. 1992| , |Cimatti et al. 1993| , 
Pi Sercgo Alighicri ct al. 1994 ) , anisotropic density 



distributions in the surrounding gas (Eales 1992) and 
inverse Compton scattering (Daly 1992). The passage 
of the radio jet still enhances star formation, but star 
formation would be not much more than in local star- 
burst galaxies, up to ~ 100M Q yr _1 . 

The possibility of a large amount of very cold molecular 
mass (temperature of a few K, as in the inner disk of 
M31, Allen & Lequeux 1993) is ruled out, because the gas 
cannot be colder than the cosmic microwave background 
which is 10 K at z ~ 2.5, so that CO lines up to J — 3— >2 



would always be relatively strong (Solomon et al. 1992b). 
Although such cold gas may show no emission from the 
higher CO transitions, the lack of CO(l-O) and CO(3-2) 
emission observed in our data indicates that even if the CO 
is cold, there must be less than ~ 10 11 M in molecular 
gas. 

It is unlikely that metallicity effects play a major role. 
The CO-H2 conversion factor a scales with CO abundance 
A (CO) approximately as a cx [A(CO)]- 1 / 4 (Radford et 
al. 1991; Solomon et al. 1992a) as long as the CO lines are 
optically thick. Abundance effects thus only begin to play 
a role when the CO lines become optically thin, which re- 
quires an abundance of several orders of magnitude below 
solar. While this situation cannot be ruled out, the pres- 
ence of large masses of dust as revealed by thermal submm 
emission in several of our targets, and the high metallici- 
ties indicated by the optical spectrum of 0211 — 122 (Van 
Ojik ct al. 1994) indicate otherwise. Simultaneously, the 
dust will shield the molecules from dissociation by the UV 



continuum, as is also the case in IRAS F10214+4724 and 
the Cloverleaf quasar (cf., Sect. 1). 

5. Summary and Conclusions 

We have searched for CO emission in a sample of 14 high 
redshift radio galaxies between z = 2 and 4.3. In none of 
these galaxies have we found consistent evidence for CO 
emission. 

For a CO-H2 conversion factor similar to the Galactic 
value, we find 2a upper limits to the molecular mass of 
less than 10 11 M . The upper limits to the molecular gas 
mass in HZRGs suggests that they are not forming stars 
at rates of ~ 10 4 M yr _1 which was previously suggested 
in some jet induced star formation scenarios to explain 
the optical-radio alignment effect. This may mean that 

(i) all powerful HZRGs have already converted their gas 
into stars, possibly through jet induced star formation, or 

(ii) jet induced star formation does not play an important 
role in the alignment effect but HZRGs may still be star 
forming at the rate of local starburst galaxies. 
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